Run

 Script to run ACA_IIR.m and CA_IIR.m



 Complex Valued Nonlinear Adaptive Filtering toolbox for MATLAB
 Supplementary to the book:
 
 "Complex Valued Nonlinear Adaptive Filters: Noncircularity, Widely Linear and Neural Models"
 by Danilo P. Mandic and Vanessa Su Lee Goh
 
 (c) Copyright Danilo P. Mandic 2009
 http://www.commsp.ee.ic.ac.uk/~mandic
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.
 
    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.
 
    You can obtain a copy of the GNU General Public License from
    http://www.gnu.org/copyleft/gpl.html or by writing to
    Free Software Foundation, Inc.,675 Mass Ave, Cambridge, MA 02139, USA.

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 % Script to run ACA_IIR.m and CA_IIR.m
0002 %
0003 %
0004 %
0005 % Complex Valued Nonlinear Adaptive Filtering toolbox for MATLAB
0006 % Supplementary to the book:
0007 %
0008 % "Complex Valued Nonlinear Adaptive Filters: Noncircularity, Widely Linear and Neural Models"
0009 % by Danilo P. Mandic and Vanessa Su Lee Goh
0010 %
0011 % (c) Copyright Danilo P. Mandic 2009
0012 % http://www.commsp.ee.ic.ac.uk/~mandic
0013 %
0014 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0015 %    This program is free software; you can redistribute it and/or modify
0016 %    it under the terms of the GNU General Public License as published by
0017 %    the Free Software Foundation; either version 2 of the License, or
0018 %    (at your option) any later version.
0019 %
0020 %    This program is distributed in the hope that it will be useful,
0021 %    but WITHOUT ANY WARRANTY; without even the implied warranty of
0022 %    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
0023 %    GNU General Public License for more details.
0024 %
0025 %    You can obtain a copy of the GNU General Public License from
0026 %    http://www.gnu.org/copyleft/gpl.html or by writing to
0027 %    Free Software Foundation, Inc.,675 Mass Ave, Cambridge, MA 02139, USA.
0028 %
0029 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0030 clear all,close all
0031 %clc
0032 % load('wind4DNord.mat')
0033 load('wind4DSud.mat')
0034 
0035 %Sampling frequency 50 Hz
0036 Freq=10;
0037 %nord=[downsample(nord(:,1),Freq) downsample(nord(:,2),Freq) downsample(nord(:,3),Freq) downsample(nord(:,4),Freq)];
0038 sud=[downsample(sud(:,1),Freq) downsample(sud(:,2),Freq) downsample(sud(:,3),Freq) downsample(sud(:,4),Freq)];
0039 
0040 beginn=1000;
0041 finish=3000;
0042 %nord=nord(beginn:finish,:);
0043 sud=sud(beginn:finish,:);
0044 
0045 %sud(:,4)= temperature, sud(:,1)= east direction, sud(:,2)= north
0046 %direction, sud(:,3)= up direction
0047 x1=[(sud(:,4)-1)*50 (sud(:,1)-1)*30 (sud(:,2)-1)*30 (sud(:,3)-1)*30];
0048 x = x1(:,2) + i*(x1(:,3));
0049 d=x;
0050 mu=1*(10^-3);
0051 M=5;
0052 
0053 %One step ahead prediction of CA-IIR ,i.e. d(n+1)
0054 [a,b,e,y] = CA_IIR(x,d,M,mu)
0055 
0056 %One step ahead prediction of Augmented CA-IIR ,i.e. d(n+1)
0057 [a,b,h,g,e,y] = ACA_IIR(x,d,M,mu)